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The exploration of novel functional carbon polymorphs is an enduring topic of scientific investigations. In 
this paper, we present simulations demonstrating metastable carbon phases as the result of pressure induced 
carbon nanotube polymerization. The configuration, bonding, electronic, and mechanical characteristics of 
carbon polymers strongly depend on the imposed hydrostatic/non-hydrostatic pressure, as well as on the 
geometry of the raw carbon nanotubes including diameter, chirality, stacking manner, and wall number. 
Especially, transition processes under hydrostatic/non-hydrostatic pressure are investigated, revealing 
unexpectedly low transition barriers and demonstrating sp^^sp^ bonding changes as well as peculiar 
oscillations of electronic property (e.g., semiconducting-^ metallic^ semiconducting transitions). These 
polymerized nanotubes show versatile and superior physical properties, such as superhardness, high tensile 
strength and ductility, and tunable electronic properties (semiconducting or metallic). 



Carbon has the ability to adopt sp-, sp^-, and sp^-hybridized states, thus forming a wide range of allotropes 
including graphite, diamond, fullerenes, carbon nanotubes (CNTs), and carbynes. The creation of new 
multifunctional carbon allotropes has been the focus of numerous theoretical and experimental explora- 
tions because of their fundamental scientific and technological importance. 

The periodic bonding motif in carbon structures can be dramatically altered under pressure, such as the well- 
known sp^^sp^ bonding change. By compressing different raw carbons such as graphite, fullerenes, CNTs, and 
vitreous carbon, diverse phases are produced including diamond, transparent superhard post-graphite phases^"^, 
one-, two-, and three-dimensional (3D) fullerene polymers^ fully sp^-bonded superhard amorphous carbon^, and 
some elusive allotropes^"^. However, completely clarifying the structures of these exciting carbon phases through 
current experimental technologies is usually unrealistic due to their limited quantity, as well as mixture of other 
phases. State-of-the-art theoretical tools for structure prediction are necessary and have yielded great success, 
such as simulated annealing^ genetic algorithm^^'^^ basin hopping^^'^^ metadynamics^^'^^, evolutionary meta- 
dynamics^^'^^, random sampling^°'^\ and particle-swarm optimization^^'^^. Most of these algorithms focus on the 
search of thermodynamically favorable structures with little consideration on the kinetics process, which is also 
important for metastable materials. The recently developed variable-cell nudged elastic band method^^ and 
transition path sampling method^^ enable further investigations into these kinetic processes. 

By considering both thermodynamics and kinetics aspects, we proposed a novel carbon allotrope, namely Cco- 
Cg^^, which can be viewed as a 3D (2,2) CNT polymer and account for the previously undetermined superhard 
carbon phase recovered from cold compression of CNT bundles^^. Further dynamical simulations demonstrated 
the easy formation of Cco-Cg by compressing armchair CNT bundles. Using transition path sampling method, we 
found that for kinetic reasons (minimal transition barrier), room-temperature compression of graphite yields M- 
carbon, which was first theoretically proposed^^'^^ and very recently unambiguously identified experimentally^^. 
These results inspired us to search and identify more new carbon allotropes under pressure with similar con- 
siderations. Numerous experimental studies on the high-pressure behavior of CNTs have been performed, most 
of which focused on the structural, electronic, and mechanical resonances of CNTs under low deformation 
ratio^^'^°. In this region, the tubes only undergo reversible polygonization, ovalization, or racetrack- shape defor- 
mations. Possible pressure-induced polymerizations of CNTs are rarely investigated except for limited cases^^"^^ 
Multi- walled CNTs (MWCNTs) become partially transparent under a non-hydrostatic pressure of 11 GPa and 
further transform into a new hard phase with complete sp^ bonding under 16 GPa^^. Under higher non-hydro- 
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Figure 1 | Top view of 3D SWCNTs polymers along the axial directions. 

1-3: 3D (4,0)-I, -II, and -III; 4-6: 3D (5,0)-I, -II, and -III; 7-9: 3D 
(6,0)-I, -II, and -III; 10-12: 3D (2,2)-I, -II, and -III; 13-15: 3D (3,3)-I, -II, 
and -III. The corresponding x, 7, and z axes of crystal structures are 
marked. 



static pressures of 24 and 35 GPa, two novel superhard carbon phases 
with respective hardness of 62-150 GPa and 58 GPa have been 
synthesized by compressing single-walled CNTs (SWCNTs) with 
Raman spectral features indicating the polymerization of CNTs^^'^^. 
Other elusive allotropes have been obtained upon simultaneous 
application of temperature and pressure^^'^^ Although these experi- 
mental results are motivating, most of the crystal structures of these 
carbon phases derived from CNTs are unsolved, and the physical 
picture underneath the pressure-induced phase transitions of nano- 
tubes remains unclear. 

In this paper, a systematic and comprehensive study about pres- 
sure-induced CNT polymerization is presented via ah initio calcula- 
tions. Significant changes in structural, electronic, and mechanical 
properties occur along with CNTs polymerization under pressure, 
yielding a variety of novel carbon allotropes strongly depending on 
the applied pressure and the nature of the CNTs precursor. The 
surprisingly low energy barrier for these transitions indicates a feas- 
ible pathway to produce new carbon allotropes with engineered 
properties. 

Results 

Polymerization under pressure. Under hydrostatic pressure, small 
zigzag (nff) and armchair {n.n) SWCNTs can easily bond covalently 
to each other by directly shortening intertube distance from 3.4 A to 
approximately 1.55 A (Fig. 1). The compressed nanotubes manifest 
diverse structures related to the CNT species and initial stacking 
configurations. For example, 3D (4,0) -I and 3D (4,0) -III carbons 
can be selectively formed due to AAA and ABA stacking of (4,0) 
SWCNTs, respectively (Fig. S2 in Supplementary Information), 
whereas the difference between 3D (4,0) -I and 3D (4,0) -II carbons 
is due to the different orientations of CNTs inside similar tetragonal 
unit cells. As a consequence, a great number of metastable carbon 
structures can be produced by compressing small SWCNTs in our 
simulations. Some previously proposed structures, such as 3D (2,2)- 
II (Bct-C/'), 3D (2,2)-III (Cco-Cg''), 3D (5,0)-IlP^ and 3D (3,3)-P^ 
have been identified along with more new structures. The transition 
pressures for small (nff) SWCNTs are usually less than 20 GPa, 
whereas those for small {n,n) SWCNTs usually exceed 20 GPa. 
This difference can be attributed to the local flattened section in 
{n,n) SWCNTs, where the n bond repulsion prevents further 
compression (see detailed discussions below). Once this repulsion 
is overcome by higher pressure, the common [2 + 2] cycloaddition 
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Figure 2 | Top view of 3D SWCNTs polymers along the axial directions. 1-4: 3D (4,0)-I, (6,0)-IV, (2,2)-II, and (8,8), recovered from (8,0), (12,0), (4,4), 
and (10,10) SWCNTs (red regions), respectively. The corresponding x, 7, and z axes of crystal structures are marked. 
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Figure 3 | (a): structural evolution from (4,4) SWCNT bundles to 3D (4,4) carbon under hydrostatic pressure of 35 GPa; (b): band gap evolution as a 
function of intertube distance d (A), where olive spheres correspond to the compression process (@ 35 GPa) and red triangles correspond to the 
decompression process (@ 0 GPa); (c): electron band structure at d= 2.304 A. 
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Figure 4 | (a-e): structural evolution from (8,0) SWCNT bundles to 3D (4,0) -I carbon under non-hydrostatic pressure of 40 GPa along x axis and z axis, 
and 36 GPa along y axis; (f): band gaps as a function of structural evolution in the compression process (olive spheres) and the decompression 
process (red triangles); (g): band structure of evolution 12 (metallic) and evolution 23 (semiconducting) in panel (f), where the magenta curves 
correspond to the bands crossing the Fermi level, and the green shadow emphasizes the band gap. 
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Figure 5 | Energy barrier curves from 3R graphite to cubic diamond, aligned (4,4) SWCNTs to 3D (4,4) carbon, and aligned (8,0) SWCNTs to 3D 
(4,0) -I carbon at hydrostatic/non-hydrostatic pressures. Non-hydrostatic A/B GPa indicates the applied pressure is A GPa along x axis and z axis, 
and B GPa along y axis. 



happens at intertube junctions, and four-membered rings appear in 
3D {n,n) SWCNT polymer structures. 

Non-hydrostatic pressure shows a much more significant effect on 
aligned large SWCNTs than on small SWCNTs. Under pressure, 
large SWCNTs have more variants and can collapse into denser 
polymers consisting of smaller SWCNTs with higher incompress- 
ibility. For example, (4,4), (8,0), and (12,0) SWCNTs are compressed 
to form 3D (2,2) (Bct-C/' and Cco-Cg"'), 3D (4,0)-I, and 3D (6,0)-IV 
carbons, respectively (See Fig. 2). Further investigations on compres- 
sion of larger nanotubes showed that more complicated structures 
can be obtained from compUcated deformations of tube cross sec- 
tions. For instance, aligned (10,10) SWCNTs under non-hydrostatic 
pressure can first transform into 3D (8,8) CNT polymer (Fig. 2), and 
then form a bond in the local graphite-like section of the structure 
under increasing pressure. The extremely complex phase behavior in 
compressed CNTs is fundamentally determined by the structural 
diversity of CNTs. It should be noted that, except for the successful 
assignment of Cco-Cs to a cold compressed CNTs^^'^^, direct com- 
parisons between our simulated CNT polymers and other synthe- 
sized phases are hard to carry out at present. Most of experimental 
studies did not provide enough structural information^ ^"^^ In addi- 
tion, the starting materials in these experiments are usually a mixture 
of different CNTs, making such a comparison even harder. A better 
control of starting SWCNT bundles, such as improved alignment 
and chirality/diameter purification, may make such comparisons 
possible in the future. 

Structural evolution. Although many experiments have suggested 
the structural transition of nanotubes under pressure, such as 
deformation, collapse, and polymerization^^'^°, the kinetic 
mechanism of transition has not been revealed up to our 
knowledge. To elaborate the kinetics process of phase transition, 
we investigated the structure evolutions of aligned (4,4) and (8,0) 
SWCNTs under hydrostatic and non -hydrostatic pressure, 
respectively (Figs. 3 and 4). For aligned (4,4) SWCNTs, adjacent 
nanotubes interact by weak van der Waals interactions at ambient 
pressure. A band gap is open at Fermi level, which is distinct from 
metalUc isolated {n,n) SWCNTs^^ but consistent with the theoretical 
prediction^" and experimental observation^^ on (n,n) SWCNT 
bundles. After applying pressure, the intertube distance in aligned 



(4,4) SWCNTs is shortened, leading to increased intertube n-n 
coupling. The (4,4) SWCNTs polymerize to form 3D (4,4)'' 
carbon, accompanied by n-^G bonding conversion in intertube 
junctions under higher pressure. Such simple structural changes 
cause unexpectedly semiconducting-^ metallic^ semiconducting 
electronic changes, as shown in Fig. 3b. When intertube distance 
decreases from standard 3.4 A to 3.054 A, enhanced intertube n-n 
coupling induces a downward shift of the green conduction band 
(Fig. 3c) near the Fermi level and closes the band gap from 0.16 eV to 
0 eV (Fig. S3). Further shortening of intertube distance continuously 
moves the green conduction lower and the red valence band higher 
across the Fermi level, resulting in a conductive behavior of the 
compressed (4,4) tubes (Figs. 3c and S3). In this "metallic" region, 
the electronic total density of states (DOS) at the Fermi level is 
determined by the green and red bands across the Fermi level and 
demonstrates a non-monotonic variation. This observation is 
reminiscent of an earlier work about the hypothetical metallic 
ThSi2- structured carbon, where the metallicity was attributed to 
short n-n contacts^\ During the final stage of compression, as 
soon as an intertube distance of 1.804 A is reached, a sudden n-^G 
intertube bonding transition occurs companied with valence electron 
localization. A band gap would re-open at the Fermi level and 
gradually increase to 0.66 eV as the length of the a-bond is 
shortened to 1.531 A. Moreover, this 3D (4,4) carbon is quen- 
chable with the band gap increasing from 0.66 eV to 0.87 eV after 
decompressing from 35 GPa to ambient pressure. 

The structural change of aligned (8,0) SWCNTs under non-hydro- 
static pressure is more complicated than that of (4,4) SWCNTs under 
hydrostatic pressure. Moreover, the corresponding electronic 
changes are more intriguing. Specifically, the aligned (8,0) 
SWCNTs undergo four major deformations (Fig. 4) under non- 
hydrostatic pressure. In the first stage, the initially circular cross- 
section of the nanotubes becomes oval-shaped (Fig. 4b), a phenom- 
enon which has been previously extensively studied^^'^'. This 
deformation is reversible and induces a semiconductor^metal 
transition in compressed (8,0) nanotubes, consistent with previous 
investigation of isolated (7,0) nanotubes^^. In the second stage, ID 
(8,0) SWCNTs start to polymerize into 2D (8,0) nanotube polymers 
(hereafter named 2D (8,0) carbon) (Fig. 4c), similar to previous 
simulations of 2D (5,0), (7,0), and (9,0) carbon'^ The 2D (8,0) 
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Figure 6 | Enthalpy (H) changes as a function of optimization steps under 
different pressures, (a) and (b): phase transitions of 3R graphite^3C 
diamond and (4,4) ^3D (4,4) carbon, respectively. Insets (a-1) and (b-1) 
show the electron density profiles corresponding to the green marks in the 
curves, while insets (a-2) and (b-2) show the electron density profiles 
corresponding to the magenta spheres in the curves. A phase transition 
only occurs after a critical pressure being applied (100 GPa for graphite and 
35 GPa for (4,4) CNT). 



carbon remains metallic because the stronger n-n coupling contri- 
butes more delocalized electrons. In the third stage, 2D (8,0) carbon 
transforms into 3D (8,0) carbon through further intertube bond 
formation (Fig. 4d). More electrons are localized in this stage, leading 
to a semiconductor with a very narrow band gap (ca. 0.02 eV). In the 
final stage, 3D (8,0) carbon is further squashed to form fully sp^- 
bonded denser 3D (4,0) -I carbon (Fig. 4e). As a result, the band 
gap rapidly increases up to 2.4 eV. This 3D (4,0) -I carbon is quench- 
able and the band gap increases to 2.8 eV after decompressing the 
pressure. This band gap oscillation is not unique, and different trans- 
ition paths would lead to different electronic property changes. 

Low transition energy barrier. Although the enthalpy of CNT 
polymer is lower than that of its counterpart CNT, a central issue 
to rationalize such phase transition is pertinent to the transition 
barrier (Fig. 5). The result of 3R graphites 3 C diamond is also 
presented in Fig. 5 for comparison. This transition can easily be 
realized by directly shortening the graphite interlayer distance and 
buckling the neighboring sheets. The calculated energy barriers 
decrease with increasing pressure. The values of 0.26 and 0.14 eV/ 
atom at 0 and 20 GPa respectively are in good agreement with 
previously calculated 0.33^^ and 0.15^^ eV/atom, indicating the 
rationality of our method (See Methods section for details). The 
pressure-induced phase transition of (4,4) SWCNT^3D (4,4) is 
also simple. The calculated transition barrier for this transition is 
notably lower than that of graphites diamond transition at the 
same pressure, implying a more feasible transition process. For 
(8,0) SWCNTs^3D (4,0)-I transition, the energy barrier at 20/18 
GPa pressure is higher than that of graphites diamond at 20 GPa. 
However, it rapidly decreases with increasing pressure and drops to 
0 eV at 40/38 GPa, whereas 90 GPa is needed to drop the barrier to 
0 eV for graphites diamond. 

In compressed graphite or (4,4) SWCNT, the highest energy-bar- 
rier points correspond to a specific structure with interlayer or inter- 
tube spacing of ca. 2.2 A. For instance, the relaxed equilibrium 
structure at 90 GPa is still graphite-like with an interlayer distance 
of about 2.2 A because this much shortened distance induces strong n 
bond repulsion to resist external pressure (Fig. 6a). At 0 K, graphite 
remains to be metastable until a critical pressure of 100 GPa is 
applied, where n-n repulsion is conquered and n-^a bonding trans- 
ition occurs accompanied by the structural change and a large 
enthalpy drop. By comparison, a critical pressure of as low as 35 
GPa is sufficient to overcome the n-n repulsion in compressed 
(4,4) SWCNTs, far less than that of graphite (Fig. 6b). These results 
clearly indicate the easy formation of new metastable carbons 
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through compressing small- diameter CNTs, compared with the case 
of graphite compression. 

Electronic and mechanical properties. Band structure calculations 
show that 3D (4,0)-II, 3D (6,0)-II, and 3D (8,8) carbons are metallic 
and other carbon polymers are semiconducting with wide direct/ 
indirect gaps from 0 eV to 3.45 eV (Table 1 and Fig. S4). CNTs 
process excellent mechanical properties, such as high stiffness and 
tensile strength, along axial direction. They are easily deformable in 
the radial direction which limits their application. Polymerization 
would greatly improve the mechanical properties in the radial 
direction without much loses along the axial one. The 3D CNT 
polymers have obviously larger Young's modulus, bulk modulus, 
and shear modulus compared with their ID counterpart (Tables 1 
and S5). We also estimated the hardness and tensile strength of 
carbon polymers based on our developed models^^"^\ The results 
indicate that 3D polymers composed of small CNTs are all 
superhard carbons, among which metallic carbons have lower 
hardness (Tables 1 and S2). Notably, our hardness model may not 
be suitable for 3D polymers composed of large CNTs with very high 
pores because the prevailing sp^-hybridized bonds are flexible and 
may not be destroyed in response to external stress. Our tensile 
strength model can be applied to all carbon polymers. The 3D 
CNT polymers have high axial tensile strength comparable to that 
of CNTs^^, and have high radial tensile strength due to sp^-bond 
bucklings along the radial directions (Tables 1, S3 and S4). The 
Poisson's ratio^^ and B/G^^ values of these carbon allotropes 
indicate that they are more ductile (less brittle) than diamond and 
lonsdaleite. An assembly of these excellent physical properties would 
make 3D CNT polymers ultimate materials for multipurpose 
applications. 

Discussion 

Carbon allotropes are notable due to the diversity of their structures 
and properties. Syntheses of these novel carbon phases depend on 
two important aspects of the transition process, namely the enthal- 
pies of the initial and final states and the transition energy barrier. 
CNTs are a suitable choice of precursors because they possess high 
enthalpy with respect to graphite. On the other hand, the transition 
barrier can be well tuned with pressure: a lower energy barrier can be 
achieved under higher pressure^^. We thus proposed an effective 
route to achieve numerous novel carbons with versatile physical 
properties, namely the pressure induced CNT polymerization. 
Specifically, the structural, electronic, and mechanical property 
changes of small SWCNTs, such as (4,0), (5,0), (6,0), (8,0), (12,0), 
(2,2), (3,3), (4,4), and (10,10) CNTs, under hydrostatic/non-hydro- 
static pressure are presented. During the transition, pecuUar conduc- 
tance oscillations in compressed CNTs are revealed in the 
calculations, an effect that can be used in nanoelectromechanical 
systems. The estimated transition barriers of small SWCNTs under 
pressure are less than that of graphite, indicating an easier synthesis 
of polymeric CNTs. The obtained 3D polymers strongly relate to the 
raw CNTs with different diameters, chiralities, stacking manners, 
and wall numbers. These novel 3D CNT polymers show excellent 
mechanical properties, such as high hardness, strength, and ductility, 
and a wide range of electronic properties ranging from metallic to 
semiconducting. Similarly exciting polymorphs may also be access- 
ible from high-pressure studies of CNTs analogs, such as boron 
nitride, silicon, germanium, and M0S2 nanotubes. As a final remark, 
we point out the possibility of achieving more exotic carbon allo- 
tropes with superior properties from compressing MWCNT bundles 
as well as mixture of distinct nanotubes, and related simulations are 
being undertaken. 

Methods 

Carbon nanotubes packed in periodic crystal lattices, such as hexagonal, 
orthorhombic, and tetragonal arrangements, with a standard intertube spacing of 3.4 



A were constructed using the Materials Studio package^*^. Structural relaxations and 
property calculations were performed based on the density functional theory (DFT) 
as implemented in the CASTEP code^*^. The Vanderbilt ultrasoft pseudopotential was 
used and the electron-electron exchange interaction was described by the local den- 
sity approximation (LDA) exchange -correlation functional of Ceperley and Alder, as 
parameterized by Perdew and Zunger (CA-PZ)^^'^*^. The plane-wave cutoff energy was 
310 eV, and a /c-point spacing {2n X 0.04 A~^) was used to generate Monkhorst-Pack 
/c-points grids for Brillouin zone sampling^^. Primitive cells were used to calculate the 
band structures and the elastic constant. Young's modulus, bulk modulus, and shear 
modulus. 

For the transition paths of 3R graphite^3C diamond and (4,4) ^3D (4,4) under 
hydrostatic pressure, the intermediate structures are determined by linearly altering 
the lattice parameters (a, b, c, a, jS, and y) from the initial structures of 3R graphite and 
(4,4) to the final structures of 3C diamond and 3D (4,4), respectively. For example, a, 
and o^are the lattice parameters of 3R graphite and 3C diamond under 30 GPa, 
respectively. After equally divide the difference Uf — ai into 50 parts, the lattice 
parameter of an intermediate structure can be defined as a„=a, + n(af — aJ/SO, 
where n running from 1 to 49 represents the n* intermediate phase. We then optimize 
the atomic positions with fixed lattice parameter at the given hydrostatic pressure 
(30 GPa). Similarly, for the transition path of (8,0) ^3D (4,0)-I under non-hydro- 
static pressure, the intermediate structures are determined by linearly altering the 
atomic fractional coordinates from (8,0) to 3D (4,0)-I. After that, the atomic frac- 
tional positions for each intermediate structure are fixed and the lattice is fully 
optimized. The corresponding enthalpies are then achieved for the optimized 
structures. The energy barrier can then be deduced from these enthalpies for specific 
applied pressure. As shown in Figure 5, thirteen typical values are selected to describe 
the trend of enthalpy and reveal the energy barrier in the transition process under 
applied pressure. 
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